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Nonlinear relativistic ionization phenomena induced by a strong laser radiation with elliptically 
■ polarization are considered. The starting point is the classical relativistic action for a free electron 

\ moving in the electromagnetic field created by a strong laser beam. The application of the relativis- 

tic action to the classical barrier-suppression ionization is briefly discussed. Further the relativistic 
version of the Landau-Dykhne formula is employed to consider the semiclassical sub-barrier ioniza- 
, tion. Simple analytical expressions have been found for: (i) the rates of the strong-field nonlinear 

/"*\ ■ ionization including relativistic initial and final state effects; (ii) the most probable value of the 

components of the photoelectron final state momentum; (iii) the most probable direction of pho- 
toelectron emission and (iv) the distribution of the photoelectron momentum near its maximum 
value. 
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I. INTRODUCTION 



Relativistic ionization phenomena induced by strong laser light have become a topic of current interest [j^-pl. In 
the nonrelativistic theory it is assumed that the electron velocity in the initial bound state as well as in the final state 
• i-h ■ is small compared with the speed of light. However, the electrons may be accelerated up to relativistic velocities in 
an intense electromagnetic field produced by modern laser devices. If the ponderomotive energy of the electron is 
^ (— i of the order of the rest energy a relativistic consideration is required. Relativistic effects in the final states become 
O • important for an infrared laser at intensities of some 10 16 Wcm~ 2 . The minimal intensity required for relativistic 
l—— effects increases by two orders of magnitude for wavelength corresponding to visible light. Ionization phenomena 
connected with relativistic final state effects have been studied for the cases of linearly and circularly polarized laser 
radiation both in the tunnel and above-barrier regimes (|]|,|7j ■ The main relativistic effects in the final state are 
\q . |Q,|U|-^] : (i) the relativistic energy distribution and (ii) the shift of the angular distribution of the emitted electrons 
£C) ' towards the direction of propagation of incident laser beam. It has been shown that a circularly polarized laser light 
produces a large amount of relativistic electrons (l]J|,|8| • On the contrary, it has been found that the ionization rate 
for relativistic electrons is very small in the case of linear polarization p . 

Relativistic effects have also to be taken into account if the binding energy E\, in the initial state is comparable 
with the electron rest energy . A relativistic formulation is necessary for the ionization of heavy atoms or singly 
or multiply charged ions from the inner K-shell. In Refs. |||| the relativistic version of the method of imaginary time 
has been employed to calculate the ionization rate for a bound system in the presence of intense static electric and 
. magnetic fields of various configurations. Analytical expressions have been found which apply to nonrelativistic bound 

systems as well as to initial states with an energy corresponding to the upper boundary of the lower continuum. 
The present paper is aimed to consider the nonlinear photoionization connected with relativistic final states velocities 
O and/or low lying initial states from a unique point of view. Further the current work extends the investigation of 
!* ■ relativistic ionization phenomena to the case of arbitrary elliptical polarization. 

Certainly this may be done in the framework of the so-called strong field approximation [Bj . In the papers of Reiss 
and of Crawford and Reiss a relativistic version of this approximation has been given for the ionization of an 

hydrogen atom with linearly and circularly polarized light. Within this approximation one calculates the transition 
amplitude between the initial Dirac state for the hydrogen atom and the final state described by the relativistic Volkov 
wave function. Coulomb corrections are neglected in the Volkov state. Therefore the final results are obtained only 
within exponential accuracy. Analytical results for the ionization rate applying to above barrier cases as well as to 
tunneling cases have been given in Refs. However, the corresponding expressions are complicated and contains 

infinite sums over all multiphoton processes. Numerical calculations are needed to present the final results. 

In contrast to the more sophisticated investigations, such as the solution of the Dirac equation or the strong field 
approximation, we are aimed to obtain simple analytical expressions. From our final formulas the explicit dependence 
of the ionization rate and of the photoelectron spectrum on the parameters, such as binding energy of the atom, field 
strength, frequency and cllipticity of the laser radiation may be understood without the need of numerical calculations. 
In this sense our approach resembles that of Popov et al. pM and of Krainov MM . 
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II. RELATIVISTIC ACTION AND CLASSICAL BARRIER-SUPRESSION IONIZATION 



Let us start with the classical relativistic action for an electron of charge e moving in the field of an electromagnetic 
plane wave with the vector potential A(t — x/c). Here and below A denotes a two-dimensional vector in the y-z plane. 
The action may be found as a solution of the Hamilton- Jacobi equation and reads 

S f {(,-^) = mAf- r --a X -- 1 + Q2 + /2 (e-eo) + ^/ [ C Ad£--4-r t A 2 dA , (1) 
I c c 2a mc z a J^ Q 2m z c l a J^ g 

where a and / = (ai, a 2 ) are constants, r = (y, z); further is £ = t — x/c, £0 is the initial value. Assuming a harmonic 
plane wave of elliptic polarization with the electric field E = F [e y cos u>£ + ge z sin u>£] we find the following expression 
for the relativistic action 

{v a / x \ 8 6 
/ 7; * + - - 7T (£ - Co) + — [01 (cos^ - coswCo) 
c 2 V c/ 2a acu 

e 2 ] 

+a 2 g (sinw£ - sinw£o)] - (g 2 - l) (sin2cj£ - sin2w£ ) ( , (2) 

aCtUl 

where the notation 8 2 = 1 + a 2 + a| + ((1 + g 2 )/2)e 2 has been introduced, the parameter e = eF/ojmc characterizes 
the strength of relativistic effects. Further the vector potential of the laser radiation has been choosen in the form 

cF cF 
A x = , Ay = sinw£, A z = g — cosw£. (3) 

W LU 

By applying the usual Hamilton- Jacobi method we take the derivative of the action Sf with respect to the constants 
ai, a 2 and a and set the result equal to new constants 8\, 82 and 83 in order to obtain the electron trajectory under 
the influence of the wave field. We obtain that the electron motion in the field and in the laboratory coordinate 
system is given by the equations (£0 = 0), 

a 2 (t + x/c) - 8 2 £ + — (ai costjg + ga 2 sinu£) + 1 - e 2 sin 2uo^ = 83, v x =c/^' ' 



cat . ce 2c . 

J/ = PH 4 cosw£, « y = 77-. jytt. i"i + esm^g} , 

a auo a(l + /(£)) 



fl . ca 2 ce . 2c 

2 = P 2 H € - 5 — sinwf, v a = {a 2 - egcosu)£} , 

a aw a(l + /(£)) 

6 2 2e 1 - a 2 

/ (0 = -o + -n («i sin w £ - 5«2 cos w£) H 5— e 2 sin 2 w£, (4) 

a z a or 

where 8\ , 82 and /?3 together with a\ , a 2 and a have to be determined from the initial conditions for position and 
velocity. Further we have introduced the notation <5 2 = 1 + a 2 + a 2 + g 2 e 2 . 

Quantum effects may be neglected, for strong enough fields, i.e., F ^S> Fb = E^/AZ (in a.u., where Eq is the 
electron energy in the initial state and Z is the effective charge of the atomic core). In this case the ionization process 
may be described by an electron trajectory given in Eqs. (|f|). In a pure classical task the constants of motion may 
be determined from the initial velocity and position of the electron at the beginning of the laser action. However, 
the initial state is given by quantum mechanics. According to a simple classical picture of ionization, the barrier 
supression ionization (BSI) [ fL2| |, the transition occurs from the bound state to that continuum state which has zero 
velocity at the time t with the phase £ of the vector potential A(£). From this condition we have to choose the 
constants as 



a = y S 2 + 2e (a\ sinw£ — ga 2 cosw£) + (1 — g 2 )e 2 c 2 sin 2 tut; , 

0!=— esino;£, a 2 = cos w£ . (5) 

The maximal ionization rate occurs at the maximum of the electric field of the laser radiation. For our choice of the 
gauge (see Eqs. @) the electric field has its maximum at the phase £ = 0. (For the sake of simplicity we neglect 
the second maximum at £ = tt). From Eqs. (^|) we conclude that the most probable final state is described by the 
constants 



2 



a = 1 , ai — , a2 = eg , 



(6) 



Two important results follow from this derivation. 

First, consider the components of the final electron drift momentum along the beam propagation, p x — c (1 — a 2 + 
al + a 2 ,) /2a, along the major axis, p y — cai, and along the small axis of the polarization ellipse, p z = ca,2, respectively. 
FromEq. (§) we see that the photoelectrons are preferably produced with the drift momentum 

eV 

Px = —^-c, p y = Q, p z = egc. (7) 

For a laser wavelength of A = 780 nm and for laser intensities of about I = (10 18 — 10 19 ) W/cm 2 , the parameter e 
is equal to: e\ = 0.65 — 2.1 for the linearly polarized wave (g = 0) and e^ = 0.46 — 1.46 for the circularly polarized 
wave (g 2 — 1). According to the classical barrier-supression ionization model the photoelectrons are emitted with a 
relativistic drift momentum ^2 ] at these laser intensities and for sufficiently large ellipticity g . On the contrary, in 
the case of linear polarization the photoelectrons have a zero drift momentum. 

Second, the angle between the electron drift momentum components along and perpendicular to the direction of 
the laser beam propagation is shifted toward the forward direction and reads 

\Pz\ \Pz\ 2 V j 

For linearly polarized laser light we obtain tan 6* = 0. For the case of circularly polarized light (where tan 8 = e/2) 
our result coincide with that of previous works |@J| . 



III. RELATIVISTIC SEMICLASSICAL APPROACH 



Consider now the process of nonlinear ionization of a strongly bound electron with a binding energy Eb comparable 
with the rest energy. Recently the ionization process in static crossed electric and magnetic fields has been considered 
HQ. The results of this paper may be applied to the ionization in laser fields only for the case of very strong fields 
e> 1. With an increasing frequency of the laser light (especially for a tentative x-ray laser )very high laser intensities 
are required to satisfy this condition. Therefore it is necessary to generalize the result of pM] to the case of nonzero 
frequencies. We consider the sub-barrier ionization. The condition to be satisfied is the opposite to the case of pure 
classical ionization, F <C Fb, in addition we have the quasiclassical condition Tilo <C Eb. No restrictions are applied 
to the parameter e. Thus we will cover both the regime of relativistic tunnel and multiphoton ionization. 

We employ the relativistic version of the Landau-Dykhne formula |^,^] . The ionization probability in quasiclassical 
approximation and with exponential accuracy reads 

TTcxexp j-|lm (S f (0; t ) + Si(t ))\ , (9) 

where Si = Eoto is the initial part of the action, 5/ is given by Eq. (Q) (or Eq. (||)). The complex initial time to has 
to be determined from the classical turning point in the complex half-plane : 

E f (to) = mc 2 \ 1 + " 2 + /2 - — ^-M(t ) + - e — A A 2 (to) \=Eo = mc 2 - E b . (10) 

la mc a Zm z c a 

Explicitely we obtain for Ao = —iujto the following relation in the case of an elliptically polarized planar wave, 

sinh 2 A - g 2 (coshA - = 7 », (11) 



where 7(a) = f]\f\ + a 2 — 2aeo, or 7 2 (a) = (1 — a) 2 ij 2 + a7 2 , with the dimensionless initial energy e — Eo/mc 2 
and the relativistic adiabatic parameter r\ = e~ x = ujmc/eF. Eq. (||) together with Eqs. @ and ( p"l| ) expresses 
the transition rate between the initial state and the final Volkov state with abitrary momentum within exponential 
accuracy. It applies for the case of sub-barrier ionization with elliptically polarized laser light. 

We are now interested in the total ionization rate. Within exponential accuracy it suffices to find the maximum of 
the transition rates between initial state and all possible final states. Aquivalently, one has to find the minimum of 
the imaginary part of the action as a function of the final state momentum. The minimization of the imaginary part 
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of the action with respect to the components of the final state momentum leads to the following boundary conditions 
11 

(x,r)(to)=Q, Im (x,r)(t = 0) =0. (12) 
From these conditions we obtain that the most probable final state is characterized by the parameters 

Q ^ 1 + 5 L( 1+/ - 29 ^-l^! 5ini2 , ), (13) 

ai = (14) 
a 2 = (.g/ry)(sinhA /A ) . (15) 

Substituting the values Ao = —iuto and a into the final state action we obtain the probability of relativistic 
quasiclassical ionization in the field of elliptically polarized laser light. Within exponential accuracy we get 

Wtxex p(-^/(7,ff,^)) . (16) 

where 

f f% n El\ = 11 4- - , 

2 7 2 a E b 2a J \ " 3 A ) 4 7 2 a 

The magnitudes a and Ao has to be taken as the solution of Eqs. ([TT| ) and (|l^). Further 7 = y/2mEbUj/eF is the 
common adiabatic Keldysh parameter from nonrelativistic theory |3| . Equation (|l^) is the most general expression 
for the relativistic ionization rate in the quasiclassical regime and for field strength smaller than the above-barrier 
threshold. It describes both the tunnel as well as the multiphoton ionization. It is the relativistic generalization of 
previous nonrelativistic results ]T^ |. 



,, P x A , 1 + ff 2 , mc 2 {l-a) 2 \ ( 2 2 tanhAo\ sinh2Ao 

fh, 9,E b )= 1 + - + — A - 1 - g + 2g — — - — . (17) 



A. Relativistic tunnel ionization 



Consider now some limiting cases. In the limit of tunnel ionization 77 <C 1 we reproduce the static result of Refs. 
and obtain the first frequency correction 



^tunnel ^ exp 



-77 (l- 5 /3) + 0(r? ), 



ao 



5a 



(18) 



where J 7 ^ = m 2 c 3 /eh = 1.32 • 10 16 V/cm is the Schwinger field of quantum electrodynamics and o?o = (eo + 
yfel + 8)/4. In the nonrelativistic regime, = Eb/mc 2 <C 1, the parameter ao = 1 — £f,/3 + e 2 /27 and the probability 
of nonrelativistic tunnel ionization including the first relativistic and frequency corrections reads 



W 



tunnel 



oc exp ■ 
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(l-5 2 /3) 



E b 



\2vnc 2 



13 
30 



7 2 (l-5 2 /3) 



(19) 



Here the first two terms in the brackets describe the familiar nonrelativistic ionization rate including the first frequency 
correction, the next two terms are the first relativistic corrections. It follows from Eq. ( |l8| ) that the account of 
relativistic effects increases the ionization rate in comparison with the nonrelativistic rate. However, even for binding 
energies of the order of the electron rest energy the relativistic correction in the exponent is quite small. In the 
"vacuum" limit Eq. ( |l9| ) results into W cx exp{— 9Fs/2F (l — 9/40t7 2 (1 — g 2 /3))}. We find a maximal deviation of 
about 18% in the argument of the exponential from the nonrelativistic formula. Here the "vacuum" limit shall not 
be confused with the pair creation from the vacuum. It is known that there are no nonlinear vacuum phenomena for 
a plane wave Jl5[ ] . In contrast to that we deal here with the ionization of an atom being in rest in the laboratory 
system of coordinates. Nevertheless, the "vacuum" limit should be considered only as the limiting result of the 
present semiclassical approach where the effects of pair production have been neglected. The polarization of the 
vacuum becomes important if the binding energy of the atom exceeds the electron rest energy. At the binding energy 
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Ef, = 2mc 2 the single particle picture employed in this paper breaks down ultimately. The electron energy is decreased 
up to the upper limit for the energy of a free positron, and the threshold energy for the production of an electron- 
positron pair becomes zero. On the contrary, for a weak relativistic initial state El « 1 we expect only a small 
influence of pair production effects on the ionization process. An appropriate consideration of vacuum polariztion 
effects can be given only in the framework of quantum electrodynamics. However, this is beyond the scope of the 
present paper. 

B. Relativistic multiphoton ionization 



Consider now the multiphoton limit rj ^> 1. In this case the parameters Ao = In (2j/^/l — g 2 ) (or Ao = lnjy/2 In 7 
for g = ±1) and a = 1 — Eb /2\q and the ionization probability in the relativistic multiphoton limit reads 

^nuiti-ph ^ exp f-^/( 7 » 1, g, E b )j , (20) 

/(7»l,fl,gt)=ln 21 -I Eb - g^±l, (21) 

v/T^2 2 8mc 2 hi27/ v ^^2 

1 E b 



/(7>l ) g, £fc )=ln27v / 21n^--- - " = = , g = ±l. (22) 

2 8mc z In 2jy/2 In 7 

Again the first two terms in the function /( 7 ^> 1, <?,£{,) reflect the nonrelativistic result |ll|, the relativistic effects 
which lead to an enhancement of the ionization probability are condensed in the third term. 

It has been shown that there is an enhancement of ionization rate in the relativistic theory for both large and small 
rj. This should be compared with the results found by Crawford and Reiss. In their numerical calculations they also 
found an enhancement of relativistic ionization rate for a circularly polarized field and for ij > 1, but for 77 <C 1 
their results suggest a strong reduction of the ionization probability Q . For the case of linearly polarized light the 
ionization rate is found to be reduced by relativistic effects 0. However, Crawford and Reiss studied the above-barrier 
ionization of hydrogen atom within the strong-field approximation. In contrast to that we have investigated the sub- 
barrier ionization from a strongly bound electron level, which yields an enhancement of the ionization rate. This 
enhancement is connected with a smaller initial time tg. As a result the under barrier complex trajectory becomes 
shorter and the ionization rate increases in comparison with the nonrelativistic theory. Figure 1 shows the relativistic 
ionization rate Eq. ([l6]) and the nonrelativistic Keldysh formula as a function of the binding energy ej, and for the 
case of linear polarization. The figure should be considered only as an illustration of the enhancement effect. The 
frequency and intensity parameters used for the calculations are still not available for the experimentalists. 



C. The case of weak relativistic initial state 



The switch from the multiphoton to the tunnel regime with increasing field strength may be studied in the nonrel- 
ativistic limit <C 1 . Within first order of the ionization probability is found to be 

^wcak-rcl K exp j_^ J(7 5 £fe <K 1} j ( (23) 

where 

/(7,5,£ 6 «1) = / (0) (7,S)+E b / (1) (7,S)- 

Here 

,rrvw x ( l + 9 2 Awm ( ? 2 tanhA (0) \ sinh2A (0) . x 

represents the nonrelativistic result p3[, and A*- ) satisfies the equation: 



smh 2 A(°) - g 2 fcoshAW - = 7 2 - 
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Besides, 
r(i) 



7 



1 + .9 2 



cosh2A^ 



.4 

- (l + .9 Z + 2 5 
is the first relativistic correction, with 



„ 2 tanhA(°)\ 2 
2 ^ ^(ol-J 



tanhA(°) 
AW 



/ sinh2A(°) \ 
{ ~ 2A(°) ) 



sinh 2 A(°) 1-5 2 . ,„, fn i\ wol A 2 „ 2 tanhA(°)\ . ,„,roil 
IPS- + "W smh2A( ] J A( 5 + (/ " 9 + 25 ^(^J sinh2A(°)| , (25) 



.l = smh2A""-2,/-(coBhAW-^^^^ 



B = 1 + g 2 



„ 2 sinh 2 A(°) 1-g 2 . , „ Wn1 
2g — ttttt^ rT7 |^sinh2A (0) . 



A(°) 2 



2A(°) 



(26) 



(27) 



Equation ( |23| ) is valid in the whole 7-domain, i.e., in the multiphoton regime 7 < 1 as well as in the tunnel limit 
7 > 1. For small adiabatic parameters, i.e., 7 — * 0, it coincides with Eq. (JlSj) ; in the case of large 7 — ► 00 it transforms 
to Eq. (|20|) . We mention that Eq. (^3|) reproduces the full relativistic formula Eq. (|l^) with very high accuracy for 
E\, < rru? . 

The expression for the rate of ionization of a weak relativistic initial state essentially simplifies in the case of linear 
polarization. Then we have 



IT' 



weak— rcl 



oc exp ■ 



2E b 



/(7,0 = O,e 6 «1) 



f(l, 9 = 0, e b < 1) = arsinh7 + 



ar 



sinh7 — 7\/i + 7 2 



- £b- 



27-^/1 + 7 2 arsinh7 + arsinh 2 7 
87 4 arsinh7 



(28) 



The terms in /(7, <? = 0, eb <C 1) which do not vanish as er& — * represent the nonrelativistic quasiclassical ionization 
rate found by Keldysh p3 ; the terms proportional to e& are the hrst relativistic correction to the Keldysh formula. 



IV. RELATIVISTIC PHOTOELECTRON SPECTRUM 



Consider now the modifications of the energy spectrum induced by relativistic effects. First we will characterize the 
most probable final state of the ejected electron. The classical nonrelativistic barrier-supression ionization predicts 
a nonzero leaving velocity of the photoelectron. However, relativistic effects as well as frequency corrections modify 
this result of the classical BSI picture. In the relativistic semiclassical theory employed in this paper we may set the 
constants a\ = and a 2 = ((//^(sinhAo/Ao) according to Eqs. ( |l4| ) and (15). From Eqs. (^) we obtain then for the 
most probable emission velocity in the laboratory system of coordinates 



vx, leaving C- 



1 " 2 +$(i-^)' 



Vy, leaving ; 
leaving = 



2ac 



g / sinhAo 



- 1 



1 + a 2 + ^ (l _ sjr^o \ A„ 

where a has to be taken from the Eq. (|l3|). In the tunnel limit (77 << 1) we obtain: 

Dx, leaving C- ■ 2 -j- Oij] ) j 

i + a 

1 — «n ^1 3n 

w 2 ,ieavin g = gcr]a —— — 2 + (J{r] ). 

1 ~r O^q 



(29) 

(30) 
(31) 



(32) 
(33) 



Here and below ao = (eo + \/ £ o + 8) /4- The first term in the x-component of the leaving velocity is independent from 
both the frequency and the intensity of the laser light. It coincides with the static result of Mur et al. Q. The leading 
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term in the z-component is proportional to the frequency and inverse proportional to the electric field strength of the 
laser radiation. From Eqs. ( ^ ) and ([53]) it also follows that the rr-component of the leaving velocity vanishes in the 
nonrelativistic limit, whereas the z-component has a nonzero nonrelativistic limit. For a nonrelativistic atom, we get: 

^leaving = ^ {l + 0{v 2 / C 2 :1 2 )} , (34) 

Reaving = ^37{l + 0(t; 2 /c 2 ,7 2 )} , (35) 



where v — y/2Eb/m is the initial "atomic" velocity of the electron. In the "vacuum" limit (ao = 1/2), we have: 

3 

^leaving = + 0(ri 2 ) , (36) 
3 

^leaving = Y^9CV + 0(f] 3 ). (37) 

It follows from these equations that a strongly bound electron has a relativistic emission velocity along the direction 
of the laser beam propagation. For a nonrelativistic initial state, <C 1, the emission velocity along the beam propa- 
gation is small. Nevertheless, the mean emission velocity seems to be the most sensitive measure of the appearance of 
relativistic effects in the initial states. In Fig. 2 the x-component of the leaving velocity is plotted versus the binding 
energy of the initial state. Though we have choosen the same parameters of the laser beam as in Fig. 1 it should be 
mentioned that the dependence of the emission velocity ^-component on the laser parameters is rather weak. The 
main parameter determining the leaving velocity along the propagation of the laser beam is the binding energy of the 
atom. 

From Eqs. (|l3|)-([l5|) we also obtain the most probable value for each component of the final state drift momentum 
(which is the full kinetic momentum minus the field momentum). We put ai = Py m /mc, a-i — p z , m /mc and 

a = (-p x , m + Jm 2 c 2 + p 2 xm +pl. m + p 2 zjn )/mc and get 



mc J 


1 - , 


2a" | 
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mc— 
V 


shXo 


Ao 



2 g 2 sinh A | 

(39) 
(40) 



The leading terms in the tunnel limit (rj << 1) read 



_ mc_g_ 
2a r\ z 



Pz,m = mc- . (42) 



V 

For a nonrelativistic initial state and within the tunnel regime (7 <^ 1) we obtain 



e 2 F 2 g 2 f l 2 g 2 + l\ r _ 
p z , m = — 9 ( 1 + -£) , (44) 



cum \ 6 

where we have given the leading terms and the first frequency corrections. 

From Eqs. (p8|)-(p0|) we easily obtain the most probable angle of electron emission. Denote by 9 the angle between 
the polarization plane and the direction of the photoelectron drift motion; and by if the angle between the projection 
of the electron drift momentum onto the polarization plane and the smaller axis of the polarization ellipse. In the 
case of a nonrelativistic atom the most probable angles read 

Vx, m = eF\g\ f g 2 + 2 r 
\p z ,m\ 2mcuj V g 2 6 



tan6» m = — r = I 1 H 2 T I ' ^ m = • v 45 ) 



We conclude that relativistic effects produce a nonzero component of the mean electron drift momentum along the 
axis of beam propagation. As a result the mean angle of electron emission is shifted to the forward direction. 
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However, in the case of linear polarization the appearance of a nonzero i-component of the photoelectron drift 
momentum is connected with relativistic effects in the initial state. The latter are typically small except the case 
of ionization from inner shells of heavy atoms. Notice that for the nonrelativistic atom the most probable value for 
the drift momentum components as well as the expression for the peak of the angular distribution coincide with the 
corresponding expressions within the BSI model (see Sec. |n|) if one neglects the frequency corrections. 

Consider now the relativistic final state spectrum, i.e., the momentum distribution near the most probable final 
state drift momentum. The calculations will be restricted to the tunnel regime 7 <C 1. Assuming weak relativistic 
effects in the initial state, Eb <C 1, and putting 5p x — {p x — p x . m ) <§; mc, Sp z — (p z — p z . m ) <C mc and p y <§; mc, one 
obtains 



W p cx H/ tunnel exp 



_7 



[Spl - 2Sp x 8p z eg + 8p 2 z (l + 2e 2 5 2 + eV/4)] 



• exp 



Pirn 7 3 (1 



to (l + eV/2) 

,,4 



7 



3 TO 



4TO 3 c 2 (l + e 2 ff 2 /2)' 



(46) 



where p^ tunnGl is the total ionization rate Eq. (tL9j) in the weak relativistic tunnel regime. In Eq. ( |4q ) only the leading 
contributions in Sp x and 5p z have been given; in the p y distribution an additional relativistic term proportional to p y 
has been maintained which becomes the leading term in the case of static fields with 7 = 0. In the non-relativistic 
limit e <C 1 and p <C c we reproduce the results of Ref. For the cases of linear (g — 0) and circular (g = ±1) 

polarization our results are in agreement with recent derivations of Krainov |6|,|8| . 

The first exponent in Eq.(|46|) describes the momentum distribution in the plane perpendicular to the major axis of 
the polarization plane. In the nonrelativistic theory (e = 0) the width of the momentum distribution in p x coincide 
with the width of the p z distribution. The relativistic effects (which are measured by eg) destroy this symmetry in 
the (x,z)-plane. The distribution of p x becomes broader, the p z distribution becomes narrower. We also mention 
the appearance of a cross term proportional to the product 6p x 5p z which is absent in the nonrelativistic theory. In 
Fig. 3 the distribution of the projection of the photoelectron drift momentum on the axis of the beam propagation 
is shown. We consider electrons which are produced in the creation of Ne 8+ (Eb = 239 eV) ions by an elliptically 
polarized laser radiation with wave length A = 1.054 ^trn, field strength 2.5 x 10 10 V/cm and ellipticity g — 0.707. The 
relativistic momentum distribution is compared with the distribution of nonrelativistic theory. From the figure we 
see that the main effect is the shift of the maximum of the momentum distribution, the broadening remains small for 
the parameters we have considered. 

The first term in the second exponent of Eq. ([l6]) determines the nonrelativistic energy spectrum for the low energetic 
electrons moving along the major polarization axis, whereas the second, relativistic term becomes important for the 
high energy tail. A detailed analysis of the photoelectron spectrum will be given elsewhere jl6| . 

In conclusion, in this paper relativistic phenomena for the ionization of an atom in the presence of intense elliptically 
polarized laser light have been considered. The cases of relativistic classical above-barrier and semiclassical sub-barrier 
ionization have been investigated. Simple analytic expressions for the ionization rate and the relativistic photoelectron 
spectrum have been obtained. These expressions apply for relativistic effects in the initial state as well as in the final 
state. We have shown that relativistic initial state effects lead to a weak enhancement of the ionization rate in the 
sub-barrier regime. The mean emission velocity has been shown to be a more sensitive measure for the appearance of 
relativistic effects in the initial state. The more important relativistic final state effects may cause a sharp increase 
of the electron momentum projection along the propagation of elliptically polarized laser light. This results in a shift 
of the most probable angle of electron emission to the forward direction. 

Finally, the expressions obtained in this paper within exponential accuracy may be improved by taking into account 
the Coulomb interaction through the perturbation theory. The results of this paper may be also used in nuclear physics 
and quantum chromo dynamics. 
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FIGURE CAPTIONS 



(Figure 1) Absolute value of the logarithm of the ionization rate — In W versus the binding energy of initial level 
eb — Eb/mc 2 . The solid line shows the relativistic rate Eq.(|l6|), the dashed line is the nonrelativistic Keldysh 
formula (Eq. (23) without the relativistic correction term). The curves are shown for a frequency to — 100 and 
an intensity / = 8.5 • 10 7 (in a.u.). 

(Figure 2) The ^-component of the emission velocity v x /c versus the binding energy of initial level &b = Eb/mc 2 . 
The emission velocity in the nonrelativistic theory is zero. The curve is shown for a frequency u) = 100 and an 
intensity I = 8.5 • 10 7 (in a.u.). 

(Figure 3) Spectrum of the electron momentum projection along the beam propagation for electrons produced in 
the creation of Ne 8+ by an elliptically polarized laser radiation with wave length A = 1.054 /im, field strength 
2.5 x 10 10 V/ 'cm and ellipticity g = 0.707; the relativistic spectrum is taken from Eq. (|46[), the nonrelativistic 
one is Eq. fiq) with e = 0. 
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